program radial
	real :: x, y, theta, r
	real, parameter :: eps = 0.2
	real,parameter :: pi = 3.14159265359
	real :: N
	print *, 'enter the number of points'; read *, N
	do theta=0.0,2.0*pi,2.0*pi/N
		!print *, f(theta), theta
		x = f(theta)*cos(theta)
		y = f(theta)*sin(theta)
		print *,x, y
	end do
contains
	real function f(theta)
	real, intent(in) :: theta
	f = (1.0 - eps**2)/(1.0 - eps*cos(theta))
	end function f
end program